MethylPCA: a toolkit to control for confounders in methylome-wide association studies

Background In methylome-wide association studies (MWAS) there are many possible differences between cases and controls (e.g. related to life style, diet, and medication use) that may affect the methylome and produce false positive findings. An effective approach to control for these confounders is to first capture the major sources of variation in the methylation data and then regress out these components in the association analyses. This approach is, however, computationally very challenging due to the extremely large number of methylation sites in the human genome. Result We introduce MethylPCA that is specifically designed to control for potential confounders in studies where the number of methylation sites is extremely large. MethylPCA offers a complete and flexible data analysis including 1) an adaptive method that performs data reduction prior to PCA by empirically combining methylation data of neighboring sites, 2) an efficient algorithm that performs a principal component analysis (PCA) on the ultra high-dimensional data matrix, and 3) association tests. To accomplish this MethylPCA allows for parallel execution of tasks, uses C++ for CPU and I/O intensive calculations, and stores intermediate results to avoid computing the same statistics multiple times or keeping results in memory. Through simulations and an analysis of a real whole methylome MBD-seq study of 1,500 subjects we show that MethylPCA effectively controls for potential confounders. Conclusions MethylPCA provides users a convenient tool to perform MWAS. The software effectively handles the challenge in memory and speed to perform tasks that would be impossible to accomplish using existing software when millions of sites are interrogated with the sample sizes required for MWAS

One CNV Discordance in <span class="italic">NRXN1</span> Observed Upon Genome-wide Screening in 38 Pairs of Adult Healthy Monozygotic Twins

Monozygotic (MZ) twins stem from the same single fertilized egg and therefore share all their inherited genetic variation. This is one of the unequivocal facts on which genetic epidemiology and twin studies are based. To what extent this also implies that MZ twins share genotypes in adult tissues is not precisely established, but a common pragmatic assumption is that MZ twins are 100% genetically identical also in adult tissues. During the past decade, this view has been challenged by several reports, with observations of differences in post-zygotic copy number variations (CNVs) between members of the same MZ pair. In this study, we performed a systematic search for differences of CNVs within 38 adult MZ pairs who had been misclassified as dizygotic (DZ) twins by questionnaire-based assessment. Initial scoring by PennCNV suggested a total of 967 CNV discordances. The within-pair correlation in number of CNVs detected was strongly dependent on confidence score filtering and reached a plateau of r = 0.8 when restricting to CNVs detected with confidence score larger than 50. The top-ranked discordances were subsequently selected for validation by quantitative polymerase chain reaction (qPCR), from which one single ~120kb deletion in NRXN1 on chromosome 2 (bp 51017111-51136802) was validated. Despite involving an exon, no sign of cognitive/mental consequences was apparent in the affected twin pair, potentially reflecting limited or lack of expression of the transcripts containing this exon in nerve/brain

Dust Devil Tracks

Dust devils that leave dark- or light-toned tracks are common on Mars and they can also be found on the Earth’s surface. Dust devil tracks (hereinafter DDTs) are ephemeral surface features with mostly sub-annual lifetimes. Regarding their size, DDT widths can range between ∼1 m and ∼1 km, depending on the diameter of dust devil that created the track, and DDT lengths range from a few tens of meters to several kilometers, limited by the duration and horizontal ground speed of dust devils. DDTs can be classified into three main types based on their morphology and albedo in contrast to their surroundings; all are found on both planets: (a) dark continuous DDTs, (b) dark cycloidal DDTs, and (c) bright DDTs. Dark continuous DDTs are the most common type on Mars. They are characterized by their relatively homogenous and continuous low albedo surface tracks. Based on terrestrial and martian in situ studies, these DDTs most likely form when surficial dust layers are removed to expose larger-grained substrate material (coarse sands of ≥500 μm in diameter). The exposure of larger-grained materials changes the photometric properties of the surface; hence leading to lower albedo tracks because grain size is photometrically inversely proportional to the surface reflectance. However, although not observed so far, compositional differences (i.e., color differences) might also lead to albedo contrasts when dust is removed to expose substrate materials with mineralogical differences. For dark continuous DDTs, albedo drop measurements are around 2.5 % in the wavelength range of 550–850 nm on Mars and around 0.5 % in the wavelength range from 300–1100 nm on Earth. The removal of an equivalent layer thickness around 1 μm is sufficient for the formation of visible dark continuous DDTs on Mars and Earth. The next type of DDTs, dark cycloidal DDTs, are characterized by their low albedo pattern of overlapping scallops. Terrestrial in situ studies imply that they are formed when sand-sized material that is eroded from the outer vortex area of a dust devil is redeposited in annular patterns in the central vortex region. This type of DDT can also be found in on Mars in orbital image data, and although in situ studies are lacking, terrestrial analog studies, laboratory work, and numerical modeling suggest they have the same formation mechanism as those on Earth. Finally, bright DDTs are characterized by their continuous track pattern and high albedo compared to their undisturbed surroundings. They are found on both planets, but to date they have only been analyzed in situ on Earth. Here, the destruction of aggregates of dust, silt and sand by dust devils leads to smooth surfaces in contrast to the undisturbed rough surfaces surrounding the track. The resulting change in photometric properties occurs because the smoother surfaces have a higher reflectance compared to the surrounding rough surface, leading to bright DDTs. On Mars, the destruction of surficial dust-aggregates may also lead to bright DDTs. However, higher reflective surfaces may be produced by other formation mechanisms, such as dust compaction by passing dust devils, as this may also cause changes in photometric properties. On Mars, DDTs in general are found at all elevations and on a global scale, except on the permanent polar caps. DDT maximum areal densities occur during spring and summer in both hemispheres produced by an increase in dust devil activity caused by maximum insolation. Regionally, dust devil densities vary spatially likely controlled by changes in dust cover thicknesses and substrate materials. This variability makes it difficult to infer dust devil activity from DDT frequencies. Furthermore, only a fraction of dust devils leave tracks. However, DDTs can be used as proxies for dust devil lifetimes and wind directions and speeds, and they can also be used to predict lander or rover solar panel clearing events. Overall, the high DDT frequency in many areas on Mars leads to drastic albedo changes that affect large-scale weather patterns

Mindfulness meditators show altered distributions of early and late neural activity markers of attention in a response inhibition task

Attention is vital for optimal behavioural performance in every-day life. Mindfulness meditation has been shown to enhance attention. However, the components of attention altered by meditation and the related neural activities are underexplored. In particular, the contributions of inhibitory processes and sustained attention are not well understood. To address these points, 34 meditators were compared to 28 age and gender matched controls during electroencephalography (EEG) recordings of neural activity during a Go/Nogo response inhibition task. This task generates a P3 event related potential, which is related to response inhibition processes in Nogo trials, and attention processes across both trial types. Compared with controls, meditators were more accurate at responding to Go and Nogo trials. Meditators showed a more frontally distributed P3 to both Go and Nogo trials, suggesting more frontal involvement in sustained attention rather than activity specific to response inhibition. Unexpectedly, meditators also showed increased positivity over the right parietal cortex prior to visual information reaching the occipital cortex (during the pre-C1 window). Both results were positively related to increased accuracy across both groups. The results suggest that meditators show altered engagement of neural regions related to attention, including both higher order processes generated by frontal regions, and sensory anticipation processes generated by poster regions. This activity may reflect an increased capacity to modulate a range of neural processes in order to meet task requirements. This increased capacity may underlie the improved attentional function observed in mindfulness meditators.Neil W. Bailey, Gabrielle Freedman, Kavya Raj, Caley M. Sullivan, Nigel C. Rogasch, Sung W. Chung, Kate E. Hoy, Richard Chambers, Craig Hassed, Nicholas T. Van Dam, Thomas Koenig, Paul B. Fitzgeral

Quantifying Between-Cohort and Between-Sex Genetic Heterogeneity in Major Depressive Disorder

Major depressive disorder (MDD) is clinically heterogeneous with prevalence rates twice as high in women as in men. There are many possible sources of heterogeneity in MDD most of which are not measured in a sufficiently comparable way across study samples. Here, we assess genetic heterogeneity based on two fundamental measures, between-cohort and between-sex heterogeneity. First, we used genome-wide association study (GWAS) summary statistics to investigate between-cohort genetic heterogeneity using the 29 research cohorts of the Psychiatric Genomics Consortium (PGC; N cases = 16,823, N controls = 25,632) and found that some of the cohort heterogeneity can be attributed to ascertainment differences (such as recruitment of cases from hospital vs community sources). Second, we evaluated between-sex genetic heterogeneity using GWAS summary statistics from the PGC, Kaiser Permanente GERA, UK Biobank and the Danish iPSYCH studies but did not find convincing evidence for genetic differences between the sexes. We conclude that there is no evidence that the heterogeneity between MDD data sets and between sexes reflects genetic heterogeneity. Larger sample sizes with detailed phenotypic records and genomic data remain the key to overcome heterogeneity inherent in assessment of MDD

Mental health indicators in Sweden over a 12-month period during the COVID-19 pandemic – Baseline data of the Omtanke2020 Study

Funding Information: This study was funded with grants from NordForsk (CovidMent, 105668 ), Horizon 2020 (CoMorMent, 847776 ), and the Karolinska Institutet . Funding Information: The Omtanke2020 study is supported by NordForsk (project No. 105668 ) and Karolinska Institute (Strategic Research Area in Epidemiology and Senior Researcher Award). We acknowledge The Swedish Twin Registry for access to contact information to participating twins. The Swedish Twin Registry is managed by Karolinska Institutet and receives funding through the Swedish Research Council under the grant no 2017-00641. The Funding Sources had no direct or indirect impact on the analysis and interpretation of the results. Publisher Copyright: © 2022 The AuthorsBackground: The ongoing COVID-19 pandemic has had an unprecedented impact on the lives of people globally and is expected to have profound effects on mental health. Here we aim to describe the mental health burden experienced in Sweden using baseline data of the Omtanke2020 Study. Method: We analysed self-reported, cross-sectional baseline data collected over a 12-month period (June 9, 2020–June 8, 2021) from the Omtanke2020 Study including 27,950 adults in Sweden. Participants were volunteers or actively recruited through existing cohorts and, after providing informed consent, responded to online questionnaires on socio-demographics, mental and physical health, as well as COVID-19 infection and impact. Poisson regression was fitted to assess the relative risk of demonstrating high level symptoms of depression, anxiety, and COVID-19 related distress. Result: The proportion of persons with high level of symptoms was 15.6 %, 9.5 % and 24.5 % for depression, anxiety, and COVID-19 specific post-traumatic stress disorder (PTSD), respectively. Overall, 43.4 % of the participants had significant, clinically relevant symptoms for at least one of the three mental health outcomes and 7.3 % had significant symptoms for all three outcomes. We also observed differences in the prevalence of these outcomes across strata of sex, age, recruitment type, COVID-19 status, region, and seasonality. Conclusion: While the proportion of persons with high mental health burden remains higher than the ones reported in pre-pandemic publications, our estimates are lower than previously reported levels of depression, anxiety, and PTSD during the pandemic in Sweden and elsewhere.Peer reviewe

A methylome-wide study of aging using massively parallel sequencing of the methyl-CpG-enriched genomic fraction from blood in over 700 subjects

The central importance of epigenetics to the aging process is increasingly being recognized. Here we perform a methylome-wide association study (MWAS) of aging in whole blood DNA from 718 individuals, aged 25–92 years (mean = 55). We sequenced the methyl-CpG-enriched genomic DNA fraction, averaging 67.3 million reads per subject, to obtain methylation measurements for the ∼27 million autosomal CpGs in the human genome. Following extensive quality control, we adaptively combined methylation measures for neighboring, highly-correlated CpGs into 4 344 016 CpG blocks with which we performed association testing. Eleven age-associated differentially methylated regions (DMRs) passed Bonferroni correction (P-value < 1.15 × 10−8). Top findings replicated in an independent sample set of 558 subjects using pyrosequencing of bisulfite-converted DNA (min P-value < 10−30). To examine biological themes, we selected 70 DMRs with false discovery rate of <0.1. Of these, 42 showed hypomethylation and 28 showed hypermethylation with age. Hypermethylated DMRs were more likely to overlap with CpG islands and shores. Hypomethylated DMRs were more likely to be in regions associated with polycomb/regulatory proteins (e.g. EZH2) or histone modifications H3K27ac, H3K4m1, H3K4m2, H3K4m3 and H3K9ac. Among genes implicated by the top DMRs were protocadherins, homeobox genes, MAPKs and ryanodine receptors. Several of our DMRs are at genes with potential relevance for age-related disease. This study successfully demonstrates the application of next-generation sequencing to MWAS, by interrogating a large proportion of the methylome and returning potentially novel age DMRs, in addition to replicating several loci implicated in previous studies using microarrays

Refinement of schizophrenia GWAS loci using methylome-wide association data

Recent genome-wide association studies (GWAS) have made substantial progress in identifying disease loci. The next logical step is to design functional experiments to identify disease mechanisms. This step, however, is often hampered by the large size of loci identified in GWAS that is caused by linkage disequilibrium (LD) between SNPs. In this study, we demonstrate how integrating methylome-wide association study (MWAS) results with GWAS findings can narrow down the location for a subset of the putative casual sites. We use the disease schizophrenia as an example. To handle “data analytic” variation we first combined our MWAS results with two GWAS meta-analyses (N=32,143 and 21,953), that had largely overlapping samples but different data analysis pipelines, separately. Permutation tests showed significant overlapping association signals between GWAS and MWAS findings. This significant overlap justified prioritizing loci based on the concordance principle. To further ensure that the methylation signal was not driven by chance, we successfully replicated the top three methylation findings near genes SDCCAG8, CREB1 and ATXN7 in an independent sample using targeted pyrosequencing. In contrast to the SNPs in the selected region, the methylation sites were largely uncorrelated explaining why the methylation signals implicated much smaller regions (median size 78bp). The refined loci showed considerable enrichment of genomic elements of possible functional importance and suggested specific hypotheses about schizophrenia etiology. Several hypotheses involved possible variation in transcription factor binding efficiencies

The Pediatric Cell Atlas: defining the growth phase of human development at single-cell resolution

Single-cell gene expression analyses of mammalian tissues have uncovered profound stage-specific molecular regulatory phenomena that have changed the understanding of unique cell types and signaling pathways critical for lineage determination, morphogenesis, and growth. We discuss here the case for a Pediatric Cell Atlas as part of the Human Cell Atlas consortium to provide single-cell profiles and spatial characterization of gene expression across human tissues and organs. Such data will complement adult and developmentally focused HCA projects to provide a rich cytogenomic framework for understanding not only pediatric health and disease but also environmental and genetic impacts across the human lifespan

Combined Whole Methylome and Genomewide Association Study Implicates CNTN4 in Alcohol Use

BACKGROUND: Methylome-wide association (MWAS) studies present a new way to advance the search for biological correlates for alcohol use. A challenge with methylation studies of alcohol involves the causal direction of significant methylation-alcohol associations. One way to address this issue is to combine MWAS data with genomewide association study (GWAS) data. METHODS: Here, we combined MWAS and GWAS results for alcohol use from 619 individuals. Our MWAS data were generated by next-generation sequencing of the methylated genomic DNA fraction, producing over 60 million reads per subject to interrogate methylation levels at ~27 million autosomal CpG sites in the human genome. Our GWAS included 5,571,786 single nucleotide polymorphisms (SNPs) imputed with 1000 Genomes. RESULTS: When combining the MWAS and GWAS data, our top finding was a region in an intron of CNTN4 (p = 2.55 × 10(-8) ), located between chr3: 2,555,403 and 2,555,524, encompassing SNPs rs1382874 and rs1382875. This finding was then replicated in an independent sample of 730 individuals. We used bisulfite pyrosequencing to measure methylation and found significant association with regular alcohol use in the same direction as the MWAS (p = 0.021). Rs1382874 and rs1382875 were genotyped and found to be associated in the same direction as the GWAS (p = 0.008 and p = 0.009). After integrating the MWAS and GWAS findings from the replication sample, we replicated our combined analysis finding (p = 0.0017) in CNTN4. CONCLUSIONS: Through combining methylation and SNP data, we have identified CNTN4 as a risk factor for regular alcohol use